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Abstract - We analyze the electronic transport through a model spin-1 molecule as a function of 
temperature, magnetic field and bias voltage. We consider the effect of magnetic anisotropy, which 
can be generated experimentally by stretching the molecule. In the experimentally relevant regime 
the conductance of the unstretched molecule reaches the unitary limit of the underscreened spin- 
1 Kondo effect at low temperatures. The magnetic anisotropy generates an antiferromagnetic 
coupling between the remaining spin 1/2 and a singular density of quasiparticles, producing a 
second Kondo effect and a reduced conductance. The results explain recent measurements in 
spin-1 molecules [Science 328 1370 (2010)]. 



^ Recent experimental advancements in measurement and 
b^ oontrol of molecular devices open the possibility of study- 
f_\ihg exotic electronic behavior in a controlled way and 
1 ' a llow for a detailed comparison with the predictions of 
_ . strongly correlated electron theories. In molecular junc- 
^..tion experiments, a single molecule is contacted to two 
•^j- metallic source and drain electrodes and its magnetic and 
i-H electronic properties can be controlled using a capacitively 
CN coupled gate electrode pQ , applying an external magnetic 
^. field or mechanically by stretching the molecule [2J. The 
***■ reduced size of the molecules leads to strong electron- 
ic electron interactions, which give rise to Coulomb block- 
^*— \ade effects, and strongly correlated electron phenomena 
• »as the spin-1 /2 Q] and the underscreened spin-1 Kondo 
,^ effects [2JS]. 



The recent observations of the underscreened spin-1 
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C3 Kondo effect [2H1] are particularly interesting because the 
" ground state of the system is not a Fermi liquid, as in the 
fully screened Kondo effect [5] , but a singular Fermi liquid 
composed by a diverging density of magnetic excitations 
at low energy, together with an asymptotically free spin 
1/2 6-8. In addition, in these experiments the ground 
state of the system can be modified and a quantum phase 
transition induced applying external fields. 

In this paper we model the experiments by Parks et 
al. [2] in which a spin-1 molecule, in the underscreened 
Kondo regime, was stretched while measuring the elec- 
tronic transport through it. We first show that an axial 



stretching of the molecule leads to a magnetic anisotropy 
term in the Hamiltonian which changes dramatically the 
low temperature transport properties. The magnetic 
anisotropy couples the remaining unscreened spin 1/2 and 
the local singular Fermi liquid excitations. As a con- 
sequence, the molecule shows a two stage Kondo effect 
to a Fermi liquid ground state. Stretching the molecule 
drives a Kosterlitz-Thouless quantum phase transition 
from a high-conductance singular Fermi liquid to a low- 
conductance Fermi liquid ground state. Applying an ex- 
ternal magnetic field parallel to the anisotropy axis, a 
crossing of the lowest lying molecular levels is induced and 
the conductance increases again. 

In the Co(tpy-SH)2 complex studied in Ref. [2], the Co 
atom is in the center of a nearly perfect octahedron of six 
N atoms. This splits the d levels of Co into three lower- 
energy t2 g and two higher-energy e g orbitals. To analyze 
the effect of stretching the molecule, we have considered 
a local model Hamiltonian which contains all interactions 
inside the d shell plus octahedral Ho and tetragonal Ht 
crystal fields (as described in Ref. [9]), and we have in- 
cluded the spin-orbit interaction XHso hi second-order 
perturbation theory. 

In absence of Hso, the ground state of the d 8 config- 
uration of Co 1+ is a B\g triplet with a hole in each e g 
orbital. The spin-orbit coupling mixes this state with sev- 
eral excited singlet and triplets which contain one e g and 
one t2 g hole. As a consequence, the S z = projection 
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state of the triplet \T, 0) is split from the non-zero pro- 
jection states \T, ±1). The energy difference is D = A 2 /, 
where / depends mainly on Ht and the Coulomb inte- 
grals which characterize the interactions inside the d shell 
[5]. This splitting can be described by a DS 2 term in 
the Hamiltonian. The tetragonal field Ht affects mainly 
the e g orbitals, which point towards the six ligand atoms. 
When the octahedron of these ligand atoms is elongated 
in the z direction (as in La2Cu04) it is more favorable en- 
ergetically to put the holes in the d x 2_ y 2 orbitals than in 
the d 3z 2_ r 2 1 leading to positive D (|T, 0} is favored), while 
when the octahedron is compressed in the same direction, 
the situation is the opposite, as in the Haldane system 
Y 2 BaNi0 5 , where D < [TO]. 

Taking the values of the Coulomb integrals and A which 
fit the low energy spectra of the neutral Ni atom (which 
has a d 8 electronic configuration as Co +1 ), an Ho splitting 
of 2 eV, and a splitting of the e g levels of 1 eV (favoring 
d x 2_ y 2 holes), we obtain D ~ 0.83 meV, a value which is 
consistent with the typical experimental observations. 

To study the electronic transport through the magnetic 
molecule we consider the effective Hamiltonian H = Hm + 
He + Hy , where 
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describes the two effective e g levels (a, b) of the molecule, 
which are relevant for the electronic transport, coupled 
through a Hund rule ferromagnetic exchange J < 0, with 
a stretching induced anisotropy D. We will focus on the 
parameter regime where the molecular ground state is in 
the spin S = 1 sector: \T,S Z ) with energy Et = 2e + 
J/4 + S z (DS z ~libH). 

The Hamiltonian of two non-interacting source and 
drain leads is given by He — X)k a £ a(k) c. 
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with 

a = L, R, and the coupling between the molecule and the 
leads is described by the last term in the Hamiltonian, 
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In order to model the experimental observations, we will 
consider that a single screening channel is relevant. A 
second channel would lead to a complete screening of the 
molecular spin, although in general at exponentially small 
temperatures [Sj ■ We assume that only one of the molecu- 
lar levels is coupled to the electrodes, and in what follows 
we take Vkflh = HcLb = 0. This can be done without loss 
of generality in the limit of large U, since other configura- 
tions are related by a level rotation. 

For symmetric hybridization to the electrodes, the con- 
ductance through the system is [IT] 
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Fig. 1: a) Zero-bias conductance through the molecule vs. 
temperature for different values of the magnetic anisotropy: 
D/Tfr 1 = 0.156,..., 0, with SD/Tjt 1 = 0.0156 between 
consecutive curves. Other parameters are: U = 0.25W, 
e = 0.125W, A = 0.035W, and J = -0.0051F with W the con- 
duction electron bandwidth and the unit of energy, b) Kondo 
temperature T^. Fit using Eq. ©, the fitting parameters are 
ci ~ 0.77, C2 ~ 1.07. Inset: Same as (a) with each curve scaled 
byTA-. 



where Paa( w ) i s the local electronic density of states on 
level a. Here, A = 27rpo(Vj 2 ) where po is the electronic 
density of states per spin of the electrodes at the Fermi 
level, the brackets denote the average over the Fermi sur- 
face, and f{ui) is the Fermi function. 

When the anisotropy term is positive D > the ground 
state of the system is a Fermi liquid JS1H1] and we can 
obtain an exact expression for zero temperature a-lcvel 
spectral density for spin a in the wide band limit [13 , 
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where n° is the charge in the level £ with spin a. Equa- 
tions © and (J3|) lead to the following expression for the 
zero temperature conductance J13H15J : 
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where Go = 2e 2 /h is the quantum of conductance. 

In absence of magnetic field n a t — nj a and g = 
sin [^ (n a + rib)] where ri£ = X^o- 71 -?- The zero- 
tempcrature conductance thus vanishes when the total 
number of electrons in the two levels of the molecule is 
even. From now on we will focus on charge sector where 
n a +Tib — 2, so that we expect a small conductance at low 
temperatures whenever the ground state is a Fermi liquid. 
In presence of a magnetic field H , the molecule will have 
a magnetization m = \(n^ — n\) and Eq. (jj) becomes 



g = sin (irm). 



(5) 



At finite temperatures we have calculated the zero-bias 
conductance using the numerical renormalization group 
(NRG) [16) H . Figure [2(a) shows the zero-bias conduc- 
tance as a function of temperature for several values of 



x For the NRG calculations wc kept up to 1600 states at each 
iteration and used a energy discretization parameter A = 2.5. 
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the anisotropy and H = 0. In absence of anisotropy, 
the molecular spin is partially screened by the conduction 
electrons below a Kondo temperature Tj| =1 . The develop- 
ment of the underscreened Kondo effect is associated with 
a monotonic increase in the conductance which, for D = 0, 
reaches the unitary limit at zero temperature. This is in 
stark contrast to what is expected for a Fermi liquid from 
Eq. (J3J). In fact, the ground state of the system is a sin- 
gular Fermi-liquid [7]. The conductance, for D = 0, is an 
universal function of T/T| =1 : G{T) = G /(T/T| =1 ) that 
differs from the one obtained for a fully screened S = 1/2 
Kondo effect |S]. 

When a positive anisotropy (D > 0) is turned on, 
the ground state of the isolated molecule is the S z = 
state \T, 0). In this case we expect a Fermi liquid ground 
state [T2] for the system, and a small conductance at low 
temperatures. For D <C ksT^ =1 we observe first an in- 
crease in the conductance as the temperature is lowered, 
followed by a plateau of high conductance and a reduction 
of the conductance for temperatures of the order of a char- 
acteristic temperature T£ <C D/ks, where T£ is defined 
by G(T/Tk) = 0.5. We note that the low-temperature 
region of the G(T) curves collapses into a single univer- 
sal curve when the temperature is scaled by T^{D) [see 
inset in Fig. |TJb)]. In the regime where D > ^flTf" 1 , 
the decrease in the conductance occurs at temperatures 
T ~ D/kB, where the energy gap D between the Sz = ±1 
and the lowest lying S z = dominates the physics and 
cuts-off the underscreened Kondo effect. 

From the plot of ln(T£) as a function of {T^ 1 / D) 1 ' 2 
[sec Fig. Gib)] it i s scen immediately that a good fit of T£ 
can be obtained using the formula 
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where in practice c\ , c 2 ~ 1 and depend weakly on the 
model parameters. This temperature scale can be identi- 
fied with the Kondo temperature of a second stage Kondo 
effect, induced by the magnetic anisotropy, in which the 
remaining spin 1/2 is screened. For D — > there is quan- 
tum phase transition of the Kostcrlitz-Thouless type from 
Fermi liquid to singular Fermi liquid. A similar behavior 
is obtained for the singlet-triplet quantum-phase transi- 
tion in models of magnetic impurities |17H19j and quantum 
dots [T4l[T5l[20H22] . Note, however, that the singlet state 
(ruled out by ab initio calculations) and a second screen- 
ing channel (ruled out in Ref. [5] by experiments with an 
external magnetic field at different angles with the stretch- 
ing axis) are absent in our model and the low-temperature 
screening is of a different nature. 

To analyze the effect of the anisotropy term, we will 
focus in the regime of large ferromagnetic \J\. In this 
regime the spin of the molecule is one, and we can describe 
it as a sum of two spin-1/2: Si and S2- Using S 2 = 
(Si + S2) 2 = 2 we can rewrite the anisotropy term as 

H D = DS 2 Z = 'W/A + D{S lz S 2z - S lx S 2x - S ly S 2y ), (7) 



which describes an anisotropic Kondo coupling. For 
D = 0, we can consider that one of the two spins is 
fully screened at temperatures T <§; T^ 1 and the other 
is asymptotically free [18] ■ When an anisotropy term 
D <C ksT^- =1 is turned on, the spin will be coupled 
through Eq. |[7J) to the singular Fermi liquid that results 
from the underscreened Kondo effect. D flows to strong 
coupling at low energies leading to a second stage Kondo 
effect. This situation is analogous to the one obtained with 
the same model, but with a positive exchange coupling J, 
where a two stage Kondo effect is obtained Q~3]. In that 
case the remaining spin couples antifcrromagnetically to 
a local Fermi liquid that results from the fully screened 
spin-1/2 Kondo effect and has a quasiparticle density of 
states: p%, P = 2+7*2 ; where Tk is the Kondo tempera- 
ture of the first Kondo stage. Performing a slave boson 
mean field approximation (SBMFT) it is then possible to 
obtain the Kondo temperature for the second stage Kondo 
effect using the formula fT3"] : 
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The resulting second stage Kondo temperature is given 
by To ~ TR-e _,rTK / J and has an excellent agreement with 
the numerical results [13] . In the present case, when the 
anisotropy term is present, the ground state of the sys- 
tem is a Fermi liquid and we expect the standard SBMFT 
approximation to give a good description of the low en- 
ergy physics. Following the same reasoning we can ex- 
tract T£ making the replacement: Tq — > T£, J — > D 
and Pq P —> Pn s P {e) in Eq. (JSJl, where Pq P is the density 
of quasiparticlcs resulting from the underscreened Kondo 
effect, to which the spin-1/2 is antiferromagnctically cou- 
pled. To recover the expression of Eq. ^ we need to 
assume the following density of quasiparticles which has 
an integrable divergence at low energy 



PQ P (e) 
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Figure [2] shows a map of the spectral density p aa (u;) as 
a function of energy and magnetic field. At H — we ob- 
serve a splitting of the spectral density of order D, since 
in this case D - 3fc B Tj| =1 . For D < fc B T| =1 , the split- 
ting is of order T^-. As the magnetic field is increased, the 
splitting is reduced and the conductance increases. The 
zero temperature conductance is proportional to p aa (0) 
and can be analyzed using the Fermi liquid relation. When 
a magnetic field is applied the molecule starts to polarize 
and the conductance increases. In the low field limit g 
increases quadratically with H , 
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were x 1S the spin susceptibility with x & 1/-D for D > 
k B T% =1 , and X « l/T£for D < k B T% =1 . 
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Fig. 2: Spectral density p aa {^>) (top) and zero-bias conductance 
(bottom) as a function of magnetic field for D ~ 3Tk ■ Other 
parameters as in Fig. [T] 




Fig. 3: a) Differential conductance versus bias voltage for dif- 
ferent values D/fc s T| =1 ~ 0, 0.36, 0.6, 0.96 (from top to 
bottom at V — 0). Other parameters are e — 0, J = 0.6VF 
and 2 A = 0.141 W b) Idem a) for D 
temperatures. 
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For a magnetic field such that /i^i? ~ D the \T, 1) 
and \T, 0) molecular levels can be tuned to be degenerate. 
In this case an orbital Kondo effect takes place and the 
conductance reaches the unitary limit g = 1 as expected 
for m = 1/2, and the spectral density shows a Kondo 
peak at w = 0. For very large H the molecule becomes 
fully polarized, m — > 1, and the conductance also vanishes. 

Finally, we turn our analysis to the out-of-equilibrium 
differential conductance at finite temperature. In the limit 
of large U and \J\, only 5 molecular states are relevant for 
the transport properties, the three projections of the spin 
I and the two projections of a spin doublet with the elec- 
tron localized at level b. This model was proposed for Tm 
impurities [33] and solved exactly [53]. We used the non- 
crossing approximation |I9j to calculate the conductance 
out of equilibrium. Figure [3] shows the differential conduc- 
tance as a function of bias voltage for different anisotropics 



as a function of the magnetic anisotropy at finite tempera- 
ture (left panel), and for a fixed anisotropy term as a func- 
tion of the temperature (right panel). As expected from 
the results of the zero temperature spectral density, the 
zero bias anomaly peak (ZBA) is split by the anisotropy 
term at low temperatures. At intermediate temperatures 
however it produces a reduction and a broadening of the 
ZBA. The results arc in good qualitative agreement with 
the experimental curves as a function of stretching and 
temperature [2]. 

In summary we have constructed a model to study 
the transport through a mechanically stretched magnetic 
molecule. We have shown that the stretching leads to 
a magnetic anisotropy term in the Hamiltonian which 
changes the low energy electronic properties of the molec- 
ular junction. The model reproduces qualitatively the ex- 
perimental behavior of the transport properties as a func- 
tion of temperature, bias voltage and magnetic field. This 
system allows for a detailed study of non-conventional 
electron behavior and magnetism in strongly correlated 
systems. 
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